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High-T c superconductors with Cu02 layers, manganites Lai-^Sr^MnOa, and cobaltites 
LaCo03 present several mysteries in their physical properties. Most of them are believed 
to come from the strongly-correlated nature of these materials. From the theoretical 
viewpoint, there are many hidden rocks in making the consistent description of the 
band structure and low-energy physics starting from the Fermi-liquid approach. Here 
we discuss the alternative method — multielectron approach to the electronic structure 
calculations for the Mott insulators - called LDA+GTB (local density approximation + 
generalized tight-binding) method. Its origin is a straightforward generalization of the 
Hubbard perturbation theory in the atomic limit and the multiband p — d Hamiltonian 
with the parameters calculated within LDA. We briefly discuss the method and focus 
on its applications to cuprates, manganites, and cobaltites. 

Keywords: LDA+GTB method; Strongly-correlated systems; Band structure. 



1. Introduction 

Fermi-liquid approach is one of the most attractive and useful in the modern con- 
densed matter physics. Its cornerstone is a quasiparticle defined as a renormalized 
electron - the one surrounded by the cloud of excitations. Such concept allows to use 
a number of the well-developed theoretical tools like diagram technique, path inte- 
grals, etc. The strong point of the Fermi-liquid approach is generality but in many 
cases it becomes a weak point. When it comes to real material-specific questions, 
the starting point for the study should depend on the particular material's charac- 
teristics. Huge leap in this direction was made in mid 60th by Hohenberg, Kohn, 
and ShauMJ^ 1 who formulated a density functional theory (DFT) . Because its start- 
ing point is the Shrodinger equation for the particular atomic arrangement, orbital, 
and spin configurations, this theory is often referred to as the 'first principles' or the 
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l ab initio' calculations. Augmented with the local density approximation (LDA) or 
the generalized gradient approximation (GGA) for the initially unknown quantity, 
exchange-correlation energy, DFT provides quantitative description of the ground 
state energy and the band structure of various atoms, molecules, and crystalline 
solids, especially containing s and p atoms (see e.gP). 

There are many cases, however, when the basic concept of the Fermi-liquid 
brakes and the system becomes a non-Fermi-liquid. Examples include temperature- 
dependent metal-insulator transitions in cobaltites and a pseudogap state in 
cuprates. In this case, one-electron approaches like LDA and GGA give qualita- 
tively wrong results. In particular, LDA failed to describe transition metal oxides 
with partially filled 3d-orbitals. The most pronounced failure is that LDA predicts 
La2Cu04 to be a metal whereas experimentally it is an insulator. The root of the 
problem is the unscreened on-site Coulomb interaction (Hubbard repulsion)^. In a 
single-band system on the mean-field level if Hubbard repulsion U is larger than 
the bandwidth W , it splits the band into two Hubbard subbands with a gap cx U . 
Spectral weight of a quasiparticle is redistributed between these subbands. At a 
half-filling, the Fermi level is inside the gap and the system is an insulator. In a 
multiorbital system, along with the Hubbard repulsion other local interactions like 
the Hund's exchange Jh and the interorbital Hubbard repulsion U' are present and 
provide a rich set of physical properties. Opening of the Hubbard gap and moreover 
the major role played by the local interactions near the half-filling are beyond the 
scope of LDA and GGA. 

There are several extensions to LDA which includes or simulates the effects of 
the on-site interactions. One of them is LDA+U^ and another one is SIC-LSDA 
(self-interaction-corrected local spin density approximation^. Both methods con- 
sider local interactions in the Hartree-Fock sense and result in the antiferromagnetic 
insulator as the ground state for La2Cu04 contrary to LDA, but the origin of the 
insulating gap is incorrect. In both LDA+U and SIC-LSDA, it is formed by the lo- 
cal single-electron states splitted by the spin or orbital polarization. Therefore, the 
paramagnetic phase above the Neel temperature TV of the undoped La2Cu04 will 
be metallic in spite of strong correlation regime U 3> W. There is one more signif- 
icant drawback in these approximations, namely, they disregard the redistribution 
of the spectral weight between the Hubbard subbands. Latter effect is incorporated 
in a different approach to ab initio calculations fo r strongl y correlated systems - 
LDA+DMFT (LDA -(-dynamical mean field theory jSSEESl 

The method is based 

on the self-consistent procedure where the LDA band structure is used to calcu- 
late the electron self-energy in DMFT. DMFT utilizes the fact that in the infinite 
dimensional limit of the Hu bbard m odel, D — > oo, the self-energy is momentum in- 
dependent, E(k, u) — > The remaining frequency dependence is exact 
in D — Y oo limit and carries very important information about dynamical corre- 
lations and Mott-Hubbard transition. On the other hand, the spatial correlations 
become crucial in low-dimensional systems like layered high-T c cuprates. That is 
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why the correct band dispersion and spectral intensities for these systems cannot 
be obtained withi n LDA+DM FT. Natural extension of this method, LDA+cluster 
or cellular DMFI 1 ^ 1 ^ 1 16 1 17 1 , and SDFT (spectral density functional theory^ 
provides momentum dependent self-energy and thus allow for the non-local corre- 
lations. 

Here we are going to discuss the alternative approach - the LDA+GTB method 
- to study the Mott-Hubbard insulators. From the very beginning the GTB (gen- 
eralized tight binding) method has been suggested to extend the microscopic band 
structure calculations to take the strong electronic correlations (SEC) into account 
in the Mott-Hubbard insulators like the transition metal oxides--. Similar to con- 
ventional tight binding (TB) approach we start with a particular local electron 
states (with all multiorbital effects, symmetry and chemistry) and then by a Fourier 
transform move to the momentum space and obtain a band structure. Because of 
SEC we can not use free electron local states, our local fermion in a d-orbital system 
is a quasiparticle given by the excitations between multielectron d n and d n±1 terms 
contrary to the conventional TB. In other words, GTB is the strongly correlated ver- 
sion of the TB method. The first computer codes and successful application of GTB 
has been developed for cuprateJ^. That version used the multiband p—d model for 
La 2 CuOjZD 

with a lot of empirical parameters in the Hamiltonian. To "cook" the 
ab initio approach, the hybrid LDA+GTB method has been developed^. After- 
ward, similar ideas have been used to stu dy t he GTB band structure of manganites 
Lai-zSr^MnOj^Sl anc } cobaltites LaCoOj 2 ^. 

The LDA+GTB may be considered as the straightforward development of the 
Hubbard atomic representation approach^ to real materials like 3d metal oxideJ^. 
Indeed, the GTB is a specific version of cluster perturbation theory (CPT) in the 
Hubbard A-operators representation . 

Later we are going to describe the method and its applications to three classes 
of materials with SEC: cuprates, manganities and cobaltites. The rest of the Review 
is organized as follows. In the next Section we provide the main ideas and technical 
steps of the LDA+GTB and consider different approximations to solve the Dyson 
equation in the A-representation. In Section [3] we discuss the LDA+GTB band 
structure of La2- K Sr x Cu04. Section |4] and Section [5] are devoted to manganites 
and cobaltites, respectively. Section [6] is the conclusion. 



2. LDA+GTB method 

The very first step in deriving the GTB method is to define an "electron" in system 
with SEC. Due to the strong interactions the free electrons are so heavily renor- 
malized, so that the new objects are unlikely to be called "the electrons" . Still, we 
would like to generate the output of the theory in terms of the Green function that 
is the one-particle property, so the "electron Green function". The goal is to find 
the electron Green function or, in other words, to define the electron. According to 
the exact Lehmann representational, at T = the electron Green function can be 
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written as 

^ y W - S2m U!-U m I 

where 0+ = £ m (A + 1) - £ (A) - /x, O" = £ (A) - B m (iV — 1) — /x, /i is the 
chemical potential, and numerators are equal to 

A m (£,w) = |(0,iV|a £(T |TO,iV + l)| 2 , 

B m (k,tu) = |(m,A-l|a £ jO,A)| 2 . 

Here, \m,N) is the m-th eigenstate of the N electron system, H\m,N) = 
E m \m,N). 

Since each single pole contribution on the right-hand side of Eq. ([!} corresponds 
to some QP, we interpret the Lehmann representation in the following way: electron 
is the linear superposition of QPs with the energies (&m) f° r electron addition 
(removal) and with the spectral weights A m (B m ). At finite temperatures, the 
Lehmann representation determines the QP as the excitation between two arbitrary 
\m,N + 1) and \n,N) terms with the energy Q mn — E m — E n and a temperature 
dependent spectral weighlP3. This definition is very clear. Unfortunately in general 
case it cannot be used straightforwardly because the exact eigenstates \m,N) and 
the eigenenergies E m are unknown. The Landau Fermi-liquid QP is a specific case 
of the Eq. ([1} with only one QP close to the Fermi level. For the free electron with 
energy e , all QP energies are equal to f2+ = f2~ = £o — A*- Now we are going 
to show that the GTB method is the perturbative realization of the Lehmann 
representation. 

As any other CPT approach, the GTB method starts with the exact diagonal- 
ization (ED) of the intracell part (H c ) of the multielectron Hamiltonian and treats 
the intercell part (H cc ) by a perturbation theory. Thus we make a realization of the 
Lehmann representation inside one unit cell with all local QP energies and spectral 
weights calculated via ED. The total LDA+GTB procedure consists of the following 
stepsPZl 

Step I: LDA. Calculation of the LDA band structure, construction of Wannier 
functions with the given symmetry, and computation of the one- and two- 
electron matrix elements of the TB Hamiltonian with the local and nearest- 
neighbor Coulomb interactions. 

Step II: ED. Separation of the total Hamiltonian H into the intra- and inter-cell 

parts, H — H c + H cc , where H c represents the sum of the orthogonal unit 

cells, H c — "Y^Hj. ED of a single unit cell term, Hf, and construction of 
/ 

the Hubbard A-operators Hubbard A-operators Xj q = \p) (q\ using the 
complete orthogonal set of eigenstates {\p}} of Hf. 
Step III: Perturbation theory. Within the A-representation, local interactions 
are diagonal and all intercell hoppings and long-range Coulomb interac- 
tion terms have the bilinear form in the A-operators. Various perturbation 
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approaches known for the Hubbard model in the X-representation can 
be used. The most general one includes treatment within the generalized 
Dyson equation obtained by the diagram technique^". 

Below we discuss each step in detail. 



2.1. Step I: LDA 

LDA provides us a set of Bloch functions | X E' A £) (A is band index) and band energies 
E\(k). For example, LDA band structure calculation for La2CuC>4 and Nd2CuC>4 
was done within the TB-LMTO-ASA (linear muffin-tin orbitals using atomic sphere 
approximation in the tight-binding) me tho d^. Using the Wannier functions ( W 
Fs) formalism^ or the NMTO method^, we ob tain single electron energies £a 
and hopping integrals Ty g A ' of the TB modeP2EH 

H = E ^ - + E E r /f c /^%^ 

+ 2 E E Vfg c /ACTi C /^3 C gA'oa c ff^ ,tr *' (^) 

f,g,\,\' U\ ,<72 ,CT3 ,CT4 

where Cf\ a is the annihilation operator in the Wannier representation of the hole 
at the site / on the orbital A and with the spin a, nf\ a = ct Ao .c/x CT - Note that 
a number and a symmetry of chosen W Fs are determined by the energy window 
that we are interested in. 

The values of Coulomb parameters Vfg are obtained by LDA supercell 
calculations^. For Cu in La2CuO-4, Hubbard parameter U and Hund's exchange 
Jh are equal to 10 eV and 1 eV, respectively^. 



2.2. Step II: Exact Diagonalization 

In transition metal (Me) oxides, the unit cell may be chosen as the MeO„ (n = 
6, 5, 4) cluster and usually there is a common oxygen shared by two adjacent cells. 
All other ions provide the electroneutrality and contribute to the high energy elec- 
tronic structure. In the low energy sector, they are inactive. Before ED calculations 
we solve the problem of nonorthogonality of the oxygen molecular orbitals of ad- 
jacent cells. For the tr-bonding of the 3d metal e„ electrons and ai„, b\ g oxygen 

104) 

orbitals this problem is solved explicitly using the diagonalization in fc-space^ 1 . We 
have used the same procedure for the t2 g orbitals. Such orthogonalization results in 
the renormalization of the hopping and Coulomb matrix elements in Eq. ([3]). Later 
we will work with the renormalized parameters. After the orthogonalization, the 
Hamiltonian |3]) can be written as a sum of intracell and intercell contributions 

H = H C + H CC , H c = ~^^Hf, H cc = ^^Hf g (4) 
/ f,g 
with orthogonal states in different cells described by if/. 
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E Nl {n) 




E x {n + 1 ) 



E (n+l) 



Fig. 1. The occupied ground term Eo(n), marked by the cross, and excited terms Ei(n) of the d n 
configuration as well as electron removal [addition] d"^ 1 [d n+1 ] sectors of the Hilbert space with 
energies Ei(n — 1) [Ei(n + 1)]. Vertical wavy lines show local Bose-type excitons. Solid lines with 
arrows show electron removal (index u) and electron addition (index c) Fermi- type excitations. 
Dashed line shows the virtual Fermi- type excitation from Eo(n + 1) to E\(n). 



The ED of Hf gives us the set of eigenstates \p) = \n p ,i p ) with the energy E p . 
Here, n p is the number of electrons per unit cell and i p denotes all other quantum 
numbers like spin, orbital moment, etc. We perform the ED with all possible excited 
eigenstates, not the Lancoz procedure. 

How to determine which configurations are relevant? They are found from the 
local electroneutrality. Let's consider LaMe03 with Me being a 3d element as an ex- 
ample. The ionic valency is La 3+ Me 3+ 03~ 2 , thus the 3d cation is Me 3+ (d 4 for Mn 3+ 
and d 6 for Co 3+ ). Due to the covalency, the ground state of the unit cell is given by 
the hybridization of d n p 6 + d n+1 p 5 configurations (in the spectroscopic notations 
d n + d n+1 L, where L means a ligand hok,^. Electron addition process results in 
d n+1 subspace of the Hilbert space with mixture of d n+1 + d n+2 L configuration. 
Similarly, electron removal results in d n ~ x subspace with d"" 1 + d n L + d n+1 L? 
mixture. Thus for stoichiometric compound, the three relevant subspaces, d™ -1 , 
d n , and d n+1 , of the Hilbert space are shown in Fig. [TJ For each subspace, the ED 
provides a set of multielectron states \n,i) with energy Ei(n), i — 0, 1, 2, N n . 
Within this set of multielectron terms, the charge-neutral Bose-type excitations 
with the energy uji = Ei(n) — Eo(n) are shown by the vertical wavy lines. Electron 
addition excitations (local Fermi-type QP) have energies J7 C! i = Ei(n + 1) — Eo(n). 
Here, index "c" means that QPs form the empty conductivity band. Similarly, 
the valence band is formed by the electron removal Fermi-type QPs with energies 
£l Vt i — E(j(n) — Ei(n — 1) [a hole creation with energy Ei(n — 1) — Eo(n) is shown 
in Fig. [Tl This multielectron language have been used in the spectroscopy, see for 
example^. The proper mathematical tool to study both the local QP and their 
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intercell hoppings is given by the Hubbard X-operatorJS, 

Xf = \p){q\. (5) 



with algebra given by the multiplication rule X pq Xy = 6 qr X ps , and by the com- 

rpp 



pleteness condition £X? P = 1. The last two equations reflects the fact that X- 



operators are the projective operators. The operator X vq describes the transition 
from initial state \q) to the final state \p), X pq \q) = \p). The important prop- 
erty of X-operators is that any local operator is given by linear combination of 

X-operators. Indeed, 6/ = 1 • 6 ■ 1 = £ \p) (pI Of \q) (q\ = £ (p\6 f \q) Xf . The 

p,i p,q 
commutation rule for Hubbard operators follows from the X-operators algebra and 

it is rather awkward. Nevertheless, if n p ~ n q is odd (±1, ±3, etc.) then the X pq is 

called quasifermionic operator; if n p — n q is even (0, ±2, etc.) the X pq is a quasi- 

bosonic one^l 

We will use a simplified notation 

Xf -> Xf™ -> Xf, (6) 

where the number m enumerate the excitation and plays the role of the QP band 
index; each pair (p, q) corresponds to some vector a(p, q) that is called "root vector" 
in the diagram technique^! In this notation, a single-electron (hole) creation opera- 
tor is given by a linear combination in the X-representation, Cf\ a = £ 7a ct (m) XT', 

rn 

l\a(p,q) = (P\cf\a \q). 

In the X-representation, the intracell part of the Hamiltonian is diagonal, 

H c = Y,(E P -n p p)Xf. (7) 

f,P 

The intercell hopping is given by 

fjtg m,m' 

where the matrix elements are 

4 W = J- T£'~llArn)7x'„{m'). (9) 

ct,A,A' 

All intraatomic d — d Coulomb interactions are included in H c and since H c is 
diagonal in X-operators, they treated exactly. The dominant part of the p — d and 
p — p Coulomb interactions is also included in H c and gives contribution to energies 
E p , while a small part of it (~ 10%) provides the intercell Coulomb interaction that 

is also bilinear in the X-operators, Hg ouL = £ £ Vf q ' p ' q ' XfX p ' q ' . 

f¥=gp,q,p' ,q' 
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2.3. Step III: Perturbation Theory 

The characteristic local energy scale is given by the effective Hubbard parameter 
U e ff = Eo(n + 1) + Eo(n — 1) — 2Eq(ji) that is given by the difference of the 
initial d n + d n and the excited d n_1 + d n+1 configurations^! The same energy 
can be obtained as the local gap between the conductivity and valence bands, 
U e ff = fi c .o — Q v ,o- Depending on the ratio of the bare Hubbard U and the charge 
excitation energy A pc [ = e p — Sd, the U e ff may represent the Mott- Hubbard gap for 
U < A p d or the charge transfer (CT) gap Ect for U > A p J33. The intercell hopping 
and interaction result in the dispersion and decrease the energy gap, E g < U e tf. The 
intercell hoppings, Eq. ©, and the non-local Coulomb interactions can be treated 
by a perturbation theory. We'd like to emphasize that in the X-representation the 
perturbation, Eq. ©, has exactly the same structure as the hopping Hamiltonian 
in the conventional Hubbard model. That is why the accumulated experience of the 
Hubbard model study in the X-representation can be used here. 

Single-electron Green function for a particle with momenta k, energy E, spin a, 
and orbital indices A and A', (E) = ( ('ckAu c^av / ) > ^ s gi verL by a linear com- 



bination of the Hubbard operator's Green functions D^™{E) = ( (X k ™ 



G£(E) = lxAmhl, a (m')DZT'(E). (10) 

mm' 

The question is how to determine X-operator's Green function. One of the regu- 
lar ways to treat perturbations is the diagram technique. For X-operators with their 
awkward algebra there is no standard Wick's theorem. Nevertheless, the general- 



ized Wick's theorem has been proved quite some time agcPSI. The first convenient 
version of the diagram technique has been formulated in Ref. The general r ules 
of the diagram technique for the X-operators are described in details in Ref. 
where the generalized Dyson equation has been obtained 

DkAE) = {(E- il m ) «W - h a (E)t{k) - t^E))^ P- ka {E). (11) 

Here Q, m = E p (n + 1) — E q (n) is the m-th QP local energy. Besides the self-energy 
matrix ^^(E), the unconventional term Pkv{E) called "the strength operator" 
appears due to the X-operators algebra. Similar term had been known for the spin- 
operator diagram technique^!. This term determines the QP oscillation strength 
(spectral weight) as well as the renormalized bandwidth. The Hubbard I solution^ 
is obtained by setting = and (Pk ia -(E)) = F rn S rnm i, where F m = 

\ ' mm' 

(X pp ) + (X qq ). We call F m "the occupation factor"; it provides a non-zero spectral 
weight for the QP excitation between at least partially filling eigenstates and gives 
zero spectral weight for excitations between empty states. The dispersion equation 
for the QP band structure of the Hubbard fermions in this case is given by 



det 



6 mn (E-n m )/F(m)-t mn (k) =0. (12) 
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The dispersion equation, Eq. (|12[) . is similar to the conventional TB dispersion 
equation, but instead of a single electron local energy e\ we have a local QP energy 
Q m . That is why we call this approach the "generalized TB method". 

The occupation numbers (X pp ) are calculated seli-consistently via the chemical 
potential equation, n e = = ^2 n /x 1 ^' 1 ' n ' l \. The change of the concentra- 

tion n e redistributes the occupation numbers and due to the occupation factors, 
F m , changes the QP band structure. 

In the GTB method, the intracell Green function can be found exactly 

G£'(E) = ]T l-fUm'tfS^DZiE), (13) 

m 

where D™ a (E) — F m j (E — TO + iS). Comparing this exact local Green function 
with the Lehmann representation in Eq. (TTJ, we can say that the electron here 
is a linear combination of local (Hubbard) fermions with QP energy fi m and a 
spectral weight |7A CT (m)| F m . It is exactly the same language as in the Lehmann 
representation. The difference is that it is realized locally and both QP energy and 
spectral weight are calculated explicitly. 

It should be stressed that the LDA+GTB bands are not the single electron 
conventional bands. There is no any single particle Schrodinger equation with the 
effective potential that gives the LDA+GTB band structure. These QP are excita- 
tions between different multielectron terms. The LDA+GTB bands depend on the 
multielectron term occupation numbers through P\t„{E) and Eko-(^) that should 
be calculated via the chemical potential equation. There is no rigid band behav- 
ior from the very beginning; the band structure depends on doping, temperature, 
pressure, and external fields. 

3. Cuprates 

In this section we present the results for cuprates with one Cu02 layer in the unit 
cell: hole doped cuprate La2- a; Sr a ;Cu04. Yet there is no self-consistent treatment 
of both electronic and magnetic structures; our LDA+GTB calculation have been 
carried out for a prescribed magnetic order. Thus, we consider separately the an- 
tiferromagnetic (AFM) phase with the long-range AFM order and the spin-liquid 
phase with the short-range AFM correlations. 

3.1. Band Structure of the Undoped La^CuO^ 

The ED of the multiband p — d Hamiltonian ^ for the Cu06 cluster which includes 
the apical oxygens, results in the following local eigenstates (in the hole represen- 
tation): 1) rih = 0, the vacuum state |0) formed by d 10 p 6 orbital configuration, 2) 
rih = 1, the spin doublets |cr, A) with different orbital symmetries. The lowest one is 
bi gi \a) , and the first excited is a\ g molecular orbital. 3) = 2. A set of two-hole 
singlets and triplets, spread in the energy region of about Ud ~ 10 eV. The lowest 
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Fig. 2. The LDA+GTB band structure of the AFM La 2 CuC>4 along the principal cuts of the 
Brillouin zoneL22l 



one is the 1 A\ singlet 15) that includes the Zhang- Rice singlet among other two- hole 
singlets. The first excited triplet \TM) (M — +1,0,-1) has the 3 Bi g symmetry. 
The total number of eigenstates is about 100. 

The next practical step is the calculation of the matrix elements, (0| Cf\ a |lcr'A), 
(la'\\cf\ a \2,i), and construction of X-operators for all single electron orbitals. 
Here, A stands for Cu-d x 2_ y 2, Cu-d z 2, O-b, O-a, or 0-p z orbital. For the AFM- 
ordered La2Cu04, we use the two-sublattice (A and B) version of the Hubbard I 
solution with two occupation factors, F m> A and F m ^s- Due to the effective molecular 
field, the local b\ g spin doublet is splitted so that at T = 0: (xV') = (-X' 



Jig spin uouoiei is spnweu so uiai &i l = u: j = \^~b 

1, (x^ = = 0. The GTB band structure and the DOS in the wide 

energy region with all excited two-hole states |2, i) have been calculated in Ref. ^21. 
The empty conductivity band is formed by only one Hubbard fermion, X® ,a . It is 
separated by the CT gap U e ff = Ect ~ 2 eV from the filled valence band. The 
valence band is formed by a large number of Hubbard fermions Xj' 21 and consists 
of a set of narrow bands with the total bandwidth about 6 eV. If we are interested 
in a smaller energy window around the Ect (for example, to study ARPES), it is 
possible to simplify the calculation by neglecting the high-energy states from both 
|2, i) and |1ct'A) sets. Then the minimal realistic basis is {|0) , \a) , IS*) , \TM}}. 

The X-representation for the fermionic operators in this basis is Cfd 2 _ io — 
uXf + 2a lx Xf, c fPba = vXf + 2a lb Xf, c fPaa = -/ a (aV2Xf° - Xf 2 °), 
c fdz2 a = 7 z (°V2Xf° - Xf 2 °), c fPza = lp (aV2Xf° - X^ T2a ). Here, a = -a 
and T2a stands for T(+l) or T(— 1) depending on the value of the spin label 
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a = ±1/2. The explicit form of the TB Hamiltonian (J3)) in this basis looks like the 
two-band singlet-triplet Hubbard model: 



f L 



Ej ^ Xj a + £2S^/ S + £2T E 



TMTM 



M 



Oct 

g 



+SS \rS<T \raS 

t fg X f X g 



2atf g {XfX™ + h.c.) 



+ t 



*J 2 ™ ) (vX / + 2a lb Xf ) + h.c.} 



+ t 



TT 

fg 



(*V2Xf os 



X 



oV2Xf° - X g 



g ) 

ctT2ct 



(14) 



The hopping parameters, Eq. ©, of the effective Hubbard model are expressed 
through the microscopic ab initio parameters t pd and t pp , tP = —2t pd Uf g 2uv — 
2t pp i> fg v 2 , tfs = -2t pd u fg 2j x j b - 2t pp is fg j b 2 , tf g = -2t pd fi fg (vj x + uj b ) - 
2t pp v fg v lb , tJJ = -^t pd \ fg 2 lalz + 2t pp v !sl l - 2t' pp \ fg 2 lplai t s f J = -^t pd ^ fglz + 
ZtppXfgla - 2£pp£/ s 7j,. Here fif g , ff g , Xf g , £/ s , and Xfg ar e the coefhcients of the 
oxygen group orbitals const ruc tion, and u, v, ~{ x , j b , la, lp, and j z are the matrix 
elements 7act(p, q) (see Ref. ^ for details). 

The QP band structure of La 2 Cu0 4 is shown in Fig.EJfor the T(0, 0)-M(tt, tt)- 
X(w, 0) — T(0,0) and X(tt,0) — Y(0, 7r) cuts of the square lattice Brillouin zone. 
Zero at the energy scale is not the Fermi level but rather fixed by the condition 

£ d w 2_ y 2 = 0. 

The top of the valence band is at the M — (ir/2, 7r/2) point, while the bottom 
of the conductivity band is at the X point. The dispersion of the valence band 
determined by the hybridization of the two bands which formed by either XJ S or 
XJ™ Hubbard fermions. The hybridization between them is provided by the t ST 
hopping matrix elements in Eq. (|14j) . These are fermionic bands, but frequently in 
the literature terms "singlet band" and "triplet band" are used. These terms reflect 
the final two-hole states involved in the QP excitations. The dominant spectral 
weight in the singlet band stems from the oxygen b\ g states, while for the bottom 
of the conductivity band it is from the d x 2_ y 2 states of Cu. 



3.2. Cascade of Lifshitz Quantum Phase Transitions in 
La 2 - x Sr x CuOi 

The low-energy Hamiltonian for I^-zSr^CuC^ (LSCO) is the t — t' — t" — J* 
model obtained via exclusion of the interband (through the charge-transfer gap) 
excitations. Here J* means that besides the Heisenberg exchange term a three-site 
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correlated hopping H3 is also included, H t -j* = H t j + H 3 , where 



E 



j. vScr \r a S 1 S9 ( vera vaa vera -v^tro^ 



tj _ \ " tfm^mg ( v aS yaa V SS v aSYSSySa\ 



Here tf g = f?^ is the hopping in the UHB, J/ ff = i^ g /U e ff is the exchange in- 
teraction due to the interband (UHB <H- LHB) hopping tf g = t ^ through the 
charge-transfer gap U e ff = Eqt, hole creation operator is now a\ a — 2aX^ a and 
its algebra is different from the bare fermion's one. The spin operators are also 
easily expressed via X-operators, Sj = Xf, Sf = [x a f a - XJ^ /2. 

Here we dropped out the triplet state \TM) because the triplet itself and the 
singlet-triplet excitations do not contribute to the near-Fermi level physics. 

The ab initio derived parameters for LSCO are (in eV) t — 0.93, t' = —0.12, 
t" = 0.15, t = 0.77, t' = -0.08, t" = 0.12, J = 0.29, J' = 0.003, J" = 0.007. 

Our approach is essentially a perturbation theory with the small parameter t/U 
contrary to the usual Fermi liquid perturbation expansion in terms of U which is 
large in cuprates. We use a method of irreducible Green functions which is similar 
to the Mori-type projection technique, with the zero-order Green function given by 
the well-known Hubbard I approximation. Beyond it there are spin fluctuations. To 
provide a description of them, the self-energy was calculated in the non-crossing 
approximation by neglecting vertex renormalization that is equivalent to the self- 
consistent Born approximation (SCBA^pM Resulting electron self-energy contains 
the space-time dependent spin correlation function C (q, ui) and results in the finite 
quasiparticle lifetime, ImS(k, ui) 7^ 0. Note that at low temperatures T < 10K the 
spin dynamics is much slower than the electron one. A typical spin fluctuation time, 



10 -9 sec, is much larger than the electronic time 10~ 13 se<P2l that is why we can 
safely neglect the time dependence of the spin correlation function, C(q, ui) —> C q . 
The self-energy becomes static, S(k, ui) — > S(k), and we have ImE = 0. Note that 
S(k, ui) here is the object completely different from the one in the Fermi liquid 
approach because the here it is build by the diagrams for the A-operators, not 
the standard Fermionic annihilation-creation operators a/ CT . In the usual Fermi liq- 
uid expansion dynamical self-energy definitely plays a crucial role in the lightly 
doped cuprates. Here, our theory starts from a different limit where the lowest or- 
der approximation is represented by the Hubbard I solution. The corrections to the 
strongly-correlated mean-field approach are small because the starting point is al- 
ready a reasonable approximation for the Mott-Hubbard insulator. That is proved 
by the small effect of the frequency dependence of the self-energy in Refs. I I '' I 
Moreover the doping-dependence of the FS is determined by ReS, and it is quali- 
tatively similar in our approaches and in the approach which properly takes ImS 
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into account3233]. 

The vertex corrections to the self-energy are small far from the spin-density wave 
or the charge-density wave instabilities, that is true for the moderate doping. Our 
approximation for the self-energy is done in the framework of the mode-coupling 
approximation w hich h as been proved to be quite reliable even for systems with 
strong interactio na l ^ . As shown in the spin-polaron treatment of the t — J model, 
the vertex corrections to the non-crossing approximation are small and give only 
numerical renormalization of the model parameters^! 

Green function (( A£ s | X^ a ^^ for a hole moving on the background of short- 
range AFM order is 

G(k, w) = (1 + ,X)/2 =5 , (15) 

u-e + n- ±±^ k - + E(k) 



where 



s(k) = — 2 —Ly 

q 

1 + x 2f k i k q 



1 — x , 1 — x t 2 



^k— q Jt 



k-q 



2 U, 



eff 



2 ^ 2 U eff 

-C q + A k _ q ^| - 1 + X % K„ 
2 q k V 2 U eff q 



Here, i k and i k are the Fourier transforms of hoppings tf g and tf g , respectively. 
The self-energy is determined by static spin correlation function Co n = (SqS~) 
and kinetic correlation function K$ n = (^oo-^^/ between sites and n. These 

correlation functions and their Fourier transforms C q and Kq represent the AFM 
short-range order and the valence-bond order, respectively. In contrast to approach 
of Ref. we calculate these correlation functions self-consistently up to n = 9 
(ninth coordination sphere) together with the chemical potential /i. To get the 
spin correlation function we also obtain the spin Green function ((X™\ X™))^ 
in a spherically-symmetric spin liquid state | 48 | 49 | w j^ n (^gz^ = q and the equal 
correlation functions for each spin component, (SqS~) = 2(S§S*) — Co n . Both 
Co n and K^ n are essentially doping-dependent and Co n decrease with the dopingH^. 
While the nearest neighbor function Coi is finite for all studied x up to x = 0.4 
with a kink at x — p* = 0.24, more distant spin correlations fall down to zero for 
x > p*. 

The calculated FS twice changes its topology with dopingj-S] see Fig. |3] Small 
hole pockets around (±7r/2, ±7r/2) points are present at small doping; then they 
increase in size and touch each other in the non-symmetric points k = ±7r(l, ±0.4) 
at x c i — p pt — 0.151. Above p op t, there are two FSs around (%, it) with outer 
being a hole-like and inner being an electron-like. The electron FS collapsed at 
%c2 = P* = 0.246, and at x > p* we have only one large hole surface around (n, ir). 
Similar conclusion on the coexistence of hole and electron FS at some intermediate 
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-It 31 ' ' ' ' ' ' ' ' " k x 



Fig. 3. (Color online) Mean-field Fermi surface transitions with doping x as calculated from poles 
of Eq. d 1 5 1) - There are two topological changes: first one between x = 0.13 and 0.16, and second 
one between 0.22 and 0.25; see Ref. S3 for detailed discussion. 



doping have been also drawn rece 

nthJSOEI] 

and earlier for the spin-density wave 

sate of the Hubbard model ^21. 

From this consideration it follows that the FS topological transitions in cuprates 
are induced by doping and they are due to the non-rigid band behavior of the 
quasiparticles in systems with SEC. According to the general Lifshitz analysis^! 
for the three dimensional (3D) system, a change of topology at the energy e = e c 
cither by appearance of a new segment (like we found at p*) or by change of its 

1/2 

connectivity (like at p op t) would result in the additional DOS, SN(e) ~ (e — e c ) 1 , 
and the change in the thermodynamic potential, 5fl ~ (ef — e c ) 5 ^ 2 (the QPT of 
the 2.5-order), where Ef is the Fermi energy. However, due to the strong anisotropy 
of electronic and magnetic properties, cuprates are quasi-2D and not isotropic 3D 
systems. The electron hopping perpendicular to the Cu02 layers in a single-layer 
cuprates like LSCO is negligibly small. 

The change of the FS topology at x c \ = p opt results in the logarithmic diver- 
gence of DOS, while the emergence of the new electron-like pocket below x C 2 = p* 
results in a step in DOS (Fig. HJ. The total DOS is a sum of the singular and 
regular contributions. We would like to stress that both logarithmic and step DOS 
singularities are in perfect agreement with the general properties of the van Hole sin- 
gularities for the 2D electrons^!. Contrary to the 3D systems, the thermodynamical 
potential for the 2D electrons has a singular contribution Sft ~ (ef ~ e c ) 2 for the 
step singularity and Sfl ^ (ef ~ e c ) 2 hi \ef ~ £ c| for the logarithmic singularity^!. 
Thus QPT at x C 2 — p* is of the second order, while at x c \ = p op t the singular- 
ity is stronger. It is immediately follows that the Sommcrfeld parameter 7 in the 
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Fig. 4. Regular, singular, and total density of states N (e — e c i) near the optimal doping e c \ = 
Sp(p pt) [top] and near the pseudogap critical point e c 2 = ep(p*) [bottom], as calculated from 
the Green function USD . Dotted curve shows the logarithmic fitting. Below p* = 0.24 (s < e c 2) 
a singular step-like contribution to the total DOS appears due to the appearance of the electron 
pocket. In the inset, the doping dependence of the superconducting critical temperature T c (x) is 
shown; the optimal doping is 0.151. Note that the energy e — e c i is the energy of holes. 



electronic heat capacity 7 = C e /T has also a singular step contribution at x < p* , 
and 5j oc In (ef — £ c ) oc In |a; — x op t \ near x c \ = p op t- Similar divergence in the spe- 
cific heat was found within the dynamical cluster approximation for the Hubbard 
modeP3. 
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The coincidence of x c \ with p opt and x C 2 with p* is not occasional. We comp are 
the superconducting critical temperature dependence T c (x) in the same modei^ 
as a function of doping and observe that T c (x) has a maximum at x op t (see inset in 
Fig. SI), which indeed equals to x c i. It is not a coincidence since like in the BCS the- 
ory the maximum in T c (x) is determined by the maximum DOS, and at x c \ we have 
a logarithmic singularity. Kinetic energy, E^ n — X^On^Om reveals a remarkable 

n 

kink at x C 2 = p* due to the change in DOS. Above p* , Eki n {p) / Eki n (p*) ~ 1 +P 
that is expected for a conventional 2D metal with the hole concentration nh = 1 +P 
and Ekin ~ £f ~ "/i- The extrapolation of this law below p* reveals that actual 
Ekin is much smaller. We associate this depleti on with the pseudogap formation 
and it fits very well with the Loram-Cooper - a simple free electron gas 

with a triangular pseudogap DOS, 



at< \ ) 9,\e-£F\ > E g 



Here £J g = J(p* — p)/p* is a doping dependent pseudogap and J is the nearest- 
neighbor exchange parameter. The observed depletion of the kinetic energy together 
with the jump in DOS relates the QPT at x C 2 to the pseudogap and the coincidence 
of x C 2 and p* is not occasional. 

The energy dependence of the electron self-energy is crucial and determines 
the Mott-Hubbard transition in the Hubbard model as was convincingly demon- 
strated byjhedynamical mean-field theory (DMFTjSH. Cluster generalization of 
DMFl P^ I ^^ is necessary to study electron correlations in a two-dimensional 
Cu02 layer where the nearest neighbor spin correlations require the momentum de- 
pendent self-energy. The cellular DMFT (CDMFT) method provides fc-dependent 
self-energy and results in the phase diagrams that have features similar to the ones 
experimentally observed in cuprated ^^^^^ l Recently, the exact diagonaliza- 
tion version of CDMFT (CDMFT+E D) was used to study the electronic structure 
of the doped Mott-Hubba rd insu latoJ^^ffl The sequence of the FS transforma- 
tions with doping in Refs. I "" |u ' 1 is very similar to ours and this is an additional 
prove of the our approach's validity at least at low temperatures and away from 
the Fermi-liquid regime. 

4. Band Structure of Manganite Lai_ x Sr x Mn03 

The most unusual solutions of Eq. (fl"2j) are those with the zero spectral weight 
(F m = 0). We know the QP energy and this QP has zero number of states. We 
call such QP the "virtual QP" . To obtain the non-zero spectral weight F m ^ the 
non-zero occupation of the initial or final states for the excitation a m = (p, q) is 
required. It may be achieved by doping, pressure, finite temperature, and external 
field. The virtual states have been found in GTB calculations for I^CuOi^and 
for FeBOaP^l. Recently they have been experimentally observed in F8BO3 by the 
IR (infra-red) spectroscop 
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GRXMGXY 0.01 0.02 0.03 
k DOS, arb.units 

Fig. 5. LDA+GTB dispersion of the in-gap and valence bands (on the left) and spin majority 
density of states (on the right) for Lai-^Sr^MnOa in the ferromagnetic phase (x = 0.3). The spin 
minority DOS in the same energy region is zero. Here, G = (0,0,0), R = (tt,tt,it), X = (it, 0,0), 
M = (tt, 7r, 0), and Y = (0, it, 0) are the main symmetry points of the Brillouin zone. 



Band structure calculations for manganites is more complicated than for 
cuprates for two reasons: a) the orbital ordering double the unit cell, b) the high 
spin values for d 4 (S = 2) and for d 4±1 (S = 5/2 for d 5 and 5 = 3/2 for d 3 ) increase 
the number of states in the Hilbcrt space. The construction of the J_-operators rep- 
resentation for LaMnC>3 within the high spin d 3 , d 4 , and d 5 configurations has been 
done in Ref. For the undoped LaMn03, the LDA+GTB calculations result in 
the QP band structure with the large CT gap Eqt ~ 2 eV. Above the top of the 
occupied valence band there is a virtual state with the activation energy Ae ~ 0.4 
eV. Doping or non-stoichiometry transforms this state into the in-gap narrow band. 
The LDA+GTB calculations for hole-doped La1_a.Sra.MnO3 for x = 0.2 -f- 0.3 re- 
sulted in the half-metallic ground state with the 100% spin polarization at the 
Fermi level in the ferromagnetic phase, that is illustrated in Fig. [SJ 

Transition from the ferromagnetic to the paramagnetic phase results in the band 
narrowing oc h cos(# ) at x = that is two times stronger than the double exchange 
model provideJ^. 



5. Finite Temperature Effect on the Electronic Structure of 
LaCo0 3 

The perovskitc-oxide LaCo03 has been studied intensely for many years due its 
unique magnetic properties and related insulator-metal transition 
A gradual appearance of the paramagnetism above 50K from the diamagnetic 
ground state in called the spin-state transition. Goodenough was the first who 
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suggested that instead of the Hund's rule dictated high spin (HS) 5 = 2 the strong 
crystalline electric field results in the low spin (LS) S=0 state for d 6 configura- 
tion of the Co 3+ ion, and the energy difference is very small with the spin gap 
A s = E(HS) - E(LS) ~ 100K. The thermal population of the HS state provides 
the sharp increase of the magnetic susceptibility \ with a maximum around 100K. 
The nature of the excited spin state of Co 3+ above the singlet 1 A\ g has been un- 
der debate (see for a recent review^J). Besides the original ^Ti g HS state with 
the ^2g e g configurations^! there were many indications on the intermediate spin 
(IS) 5 = 1, 3 Ti g , state. The two stage model has b een pr oposed with the LS-IS 
transition at 100K and IS-HS transition at 550- 600K™. Recent electronic spin 
resonance (ESR}P^, X-ray absorption spectroscopy (XAS) and X-ray magnetic cir- 
cular dichroism (XMCD/-^ experiments prove that the lowest excited state is really 
the HS. Nevertheless the ^Ti g term is splitted by the spin-orbital interaction in the 
low energy triplet with effective moment J = 1, and higher energy sublevels with 
J = 2 and J = 3^1 The large difference between the spin excitation gap As and 
the charge gap given by the activation energy for electrical conductivity E a w O.leV 
at low T indicates that LaCoC>3 is not a simple band insulator^! The second shal- 
low maximum in \ near 500 -j- 600 is often related to the insulator-metal transition. 
Surprisingly for the IMT electrical con ductivity a does not seem to show any notice- 
able anomaly at this temperature^^. Moreover the discrepancy between the large 
charge gap 2E a s» 2300K and the Tjmt ~ 600K implies that the IMT cannot be 
simply argued in terms of a narrow-gap semiconductor^!. We solved this problem 
by calculating the electronic band structure in the regime of strong electron corre- 
lations. We consider electron as the linear combination of quasiparticles (QP) given 
by excitations between the different multielectron configurations obtained by ex- 
act diagonalization of the CoOg cluster. With the Hubbard operators constructed 
within the exact cluster eigenstates we can calculate the QP band structure for 
the infinite lattice. The QP spectral weight is determined by the occupation num- 
bers of the local multielectron configurations. We find that the thermal population 
of different sublevels of the 5 T2 g HS term splitted by the spin-orbital interaction 
results both in the spinstatc transition and also in some new QP excitations. Of 
particular importance is the hole creation QP from the initial d 6 HS into the d 5 
HS term, this QP appears to form the in-gap state inside the large charge-transfer 
gap Eg s» 1.5eV. The intercluster hopping transforms this local QP into the in-gap 
band that lies just under the bottom of empty conductivity band and provides the 
insulating gap 2E a s» 0.2eV at T = 100K. It bandwidth increases with T, and 
overlapping with the conductivity band at T — T\mt — 587K results in the IMT. 
Hence our approach allows to treat both the low T spin-state transition and the 
high T IMT on the same footing. 

LaCoOa as well as other strongly correlated oxides is a difficult problem for 
the ab initio band theory. The LDA calculations^^ incorrectly predict a metal 
for paramagnetic LaCoOa. Various methods have been applied to study effect of 
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(a) N c =5 (b) N e =6 (c) N e =7 



Fig. 6. The low-energy part of the Hilbert space for CoOg cluster with the electron numbers 
N e = 5, 6, 7. Terms with a given N e are the mixtures of d e , d Nc + 1 L, and d Ne + 2 L 2 configurations. 
At T = 0, only the N e = 6 low-spin term 1 Ai is occupied; the Fermi- type excitations from this 
term which form the top of the valence band (d e d 5 ) and the bottom of the conductivity band 
(a! 6 — > d 7 ) are shown by the solid lines with arrows. The dashed lines denote the in-gap excitations 
with the spectral weight increasing with temperature due to the population the HS excited d 6 
terms. 



correlations on the LaCoC>3 electronic structure: LDA+U or 

dynamical me an- field theory 2 ^. Recent variational cluster approximation (VCA) 
calculation 2 ^ based on the exact diagonalization of the CoC>6 cluster gives a rea- 
sonably accurate description of the low temperature properties: the insulating na- 
ture of the material, the photoelectron spectra, the LS-HS spin-state transition. 
The main deficiency of the VCA is the failure to reproduce the high temperature 
anomalies in the magnetic and electronic properties associated with the IMT. 

The exact diagonalization of the multielectron Hamiltonian for a finite cluster 
provides a reliable general overview of the electronic structure of the correlated 

IOQI 

materials 221 . To incorporat e the lattice effect several versions of the cluster per- 
turbation theory are 

known™ 

To calculate the band dispersion in the strongly 
correlated material one has to go beyond the local multielectron language. The nat- 
ural tool to solve this problem is given by the Hubbard A-operators X^ q = \p) (q\ 
constructed with the Co06 cluster eigenvectors \p) at site Rf. All effects of the 
strong Coulomb interaction, spin-orbit coupling, covalence and the crystal field in- 
side the Co06 cluster are included in the set of the local eigenstates E p . Here p 
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Fig. 7. Quasiparticle dispersion and density of states at two different temperatures, T = 100K (a) 
and T = 600K (b). At T = OK, LaCo03 is the charge-transfer insulator with the gap E g 1.5eV. 
At finite temperatures, the in-gap band appears below the conductivity band with the temperature 
dependent activation energy. At T = fOOK, E a « O.leV (a). At T = Timt = 587K, E a = OeV, 
and above the Tjmt the band structure is of the metallic type (b). 



denotes the following quantum numbers: the number of electrons (both 3d Co and 
p of O), spin S and pseudoorbital moment I (or the total pseudomoment J due to 
spin-orbit coupling), the irreducible representation in the crystal field. A relevant 
number of electrons is determined from the electroneutrality, for stoichiometric 
LaCoOs, n = 6. In the pure ioni c m odel the corresponding energy level scheme 
for d 6 Co 3+ has been obtained irPSI. Due to the covalence there is admixture of 
the ligand hole configurations d n+1 L and d n+2 L 2 that is very well known in the 
X-ray spectroscopy^. Contrary to spectroscopy the electronic structure calcula- 
tions require the electron addition and removal excitations. For LaCoOs it means 
the d 5 and d 7 configurations. The total low energy Hilbert space is shown in the 
Fig. [51 Here the energy level notations are the same as in the ionic model^but all 
eigenstate contains the oxygen hole admixture due to the covalence effect. The cal- 
culation of the n = 5, 6, 7 eigenvectors for CoOfiCluster with the spin-orbit coupling 
and the Coulomb interaction has been done in^l 

The electron removal spectrum determines the top of the valence band, the 



I 



I 



August 23, 2012 0:20 WSPC/INSTRUCTION FILE review2012-mplb 



Cuprates, Manganites, and Cobaltites: Multielectron Approach to the Band Structure 21 

corresponding electron quasiparticle QP are shown in the Fig.|6]by thin solid lines 
as the excitation from the 1 A\ d 6 singlet in the 2 T 2 d 5 states with J = 1/2 and 
J = 3/2. There energies are 

n vl = E (d 6 , %) - E (d\ 2 T 2 , J = 1/2) , (17) 

n V2 = E (d 6 , X A 1 ) - E (d\ 2 T 2 , J = 3/2) . (18) 

The bottom of empty conductivity band has the energy 

n c =E(d 7 , 2 E)-E(d 6 , 1 A 1 ). (19) 

All these bands have nonzero QP spectral weight. The intercluster hopping results 
in the dispersion, f2„ — > J7„(fc). The QP LDA+GTB band structure corresponds to 
the insulator with the gap E g ss 1.5eV (Fig. [JJ at T = 0. This gap value is rather 
close to the VCA gap^ and the experimental value E g « leVSQI. 

At finite temperature, the thermal excitation over the spin-gap A5 into the 
J = 1 and over the gap A5 + 2A into the J = 2 sublevels of the HS 5 T 2g state 
occurs. We take As = 140K and A = — 185K following Ref. [76] Partial occupation 
of the excited HS states results in the drastically change of the QP spectrum. 
For T — OK, excitations from the l A\ d 6 singlet in the lowest 6 Ai d 5 term were 
forbidden due to spin conservation (the corresponding matrix element 7„ =0), and 
the excitation from <i 6 , J = 1^ in |<i 5 , & A\j has nonzero matrix element (shown 
by dashed line fly 1 in Fig. [6]) but zero filling factor as the excitation between two 
empty states. For T 7^ 0, the filling factor for the fly 1 and fl V2 QP is nonzero 



and is equal to the occupation number ri\ and n 2 of the states 
d & , J — 2) correspondingly. The energies of these QP are 



d 6 , J = 1 > and 



tl vl = E (d 6 , 5 T 2g , J = lj - E (d 5 , 6 A,) , (20) 
tl V2 = E (d 6 , 5 T 2g , J = 2)-E (d 5 , 6 Ax) . (21) 



Energies of these QP appear to be slightly below the bottom of the conductivity 
band, see DOS at finite temperature in Fig. [7l Thus we have obtained that these 
temperature-induced QP states lies inside the charge-transfer gap, they are the 
in-gap states. Similar in-gap states are known to results from doping in the high 
temperature superconductors. The LaCoOs is unique because the in-gap states are 
induced by heating. The chemical potential lies in the narrow gap 2E a 0.2eV 
at T = 100K between the in-gap states and conductivity band. 

From the GTB dispersion equation (fT2"|) it is clear that the in-gap bandwidth 
is proportional to the occupation numbers n\ and n 2 of the excited HS states. 
With further temperature increase, the in-gap bands and Q V2 become wider 
and finally overlap with the conductivity band (Fig. EJ) at T = Timt = 587K. 
It should be clarified that the IMT in LaCoOs is not the thermodynamic phase 



I 
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transition, there is no any order parameter associated with the gap contrary to the 
classical IMT in V0 2 , NiS etc. 

Thus, we find that a correct definition of the electron in strongly correlated 
system directly results in the in-gap states during the spin-state transition due to 
the thermal population of the excited HS states. Close to the spin-state temperature 
region the in-gap states determine the value of the activation energy E a w O.leV. 
Further temperature increase results in large in-gap bandwidth and smaller E a , and 
finally E a = at Timt- As concerns the weak maximum in the x(T) close to the 
IMT, it may be a small Pauli-type contribution from the itinerant carriers above 
Timt- We emphasize that instead of rather large difference in temperatures of the 
spin-state transition (~ 100K) and the IMT (600K) the underlying mechanism is 
the same and is induced by the thermal population of the excited HS states. 

6. Conclusions 

We have presented the main ideas of the LDA+GTB method. Being invented to 
study the high-T c superconductivity in cuprates, LDA+GTB method appears to be 
powerful approach to systems with SEC and useful for other Mott insulators. Since 
it is a combination of the ab initio and model approaches, the method cannot go 
beyond the restriction of the model used in GTB. For example, the absence of the 
long-range Coulomb interaction which determines the Coulomb matrix elements 
in the large-wavelength limit prevents the correct description of the overdoped 
cuprates. Of course, it is the common deficiency of all Hubbard-type models. The 
modern version of the LDA+GTB cannot be used when the perturbation parameter 
t/U increases and the Mott transition is expected. Nevertheless, it works in the most 
difficult for conventional band theory region of strong electronic correlations. 

Application of the LDA+GTb method for cuprates revealed two critical points in 
the doping dependence. The first one is related to the change of the FS connectivity 
and logarithmic divergences of DOS and of electronic heat capacity parameter 7 
at the optimal doping p opt = 0.151. The logarithmic enhancement of DOS leads to 
the maximum in the doping dependence of superconducting critical temperature 
T c at the same critical point x = p op t- The second QPT is associated with the 
collapse of the electron-like FS pocket at p — > p* = 0.246 and results in the step 
singularities in DOS and in Sommerfeld parameter 7. We have found the depletion 
of the hole's kinetic energy below p* and ascribe it to the pseudogap formation at 
p < p* . Thus the two energy scales in cuprates measured by T c and T* are both 
related to the QPTs and to the changes of the cuprate's electronic structure with 
doping. The underlying physics is tightly connected with the scattering on the spin 
fluctuations. While approach starting from the conventional Feynman diagrammatic 
expansion requires frequency-dependent self-energy to get the same results (see, e.g., 
CDMFT), A-operators technique allows to catch the crucial effects of scattering 
even on the mean-filed level. 

The multiorbital extension of the same ideas results in a qualitatively correct 
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description of the main peculiarities of the electronic properties of manganites 
and cobaltites. Further study of these and other Mott-Hubbard insulators are in 
progress. 
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